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Abstract — It is known that the frame error rate of turbo codes 
on quasi-static fading channels can be accurately approximated 
using the convergence threshold of the corresponding iterative 
decoder. This paper considers quasi-static fading channels and 
demonstrates that non-iterative schemes can also be character- 
ized by a similar threshold based on which their frame error rate 
can be readily estimated. In particular, we show that this thresh- 
old is a function of the probability of successful frame detection 
in additive white Gaussian noise, normalized by the squared in- 
stantaneous signal-to-noise ratio. We apply our approach to un- 
coded binary phase shift keying, convolutional coding and turbo 
coding and demonstrate that the approximated frame error rate 
is within 0.4 dB of the simulation results. Finally, we introduce 
performance evaluation plots to explore the impact of the frame 
size on the performance of the schemes under investigation. 

I. Introduction 

The performance analysis of transmission schemes on quasi- 
static fading channels constitutes an important problem owning 
to the fact that this channel model characterizes practical set- 
tings that experience extremely slow fading conditions, such 
as fixed wireless access systems [1]. In this context, bounding 
techniques for the error rate of various transmission schemes 
have been proposed; such schemes include block codes [2], 
convolutional codes [3], [4], turbo codes [5], space-time trellis 
codes [6], [7] and serially concatenated codes [5], [8]. 

El Gamal and Hammons [9] have proposed an analytical ap- 
proximation to the frame error rate (FER) of iterative schemes, 
such as turbo codes, on quasi-static fading channels; this tight 
approximation is made possible owning to a simple character- 
ization of an iterative decoder. In this paper, we will demon- 
strate that non-iterative schemes, both uncoded and coded, can 
also be characterized in a similar manner, hence their FER 
performance can be estimated using the same approximation 
technique. 

The rest of the paper is organized as follows. The quasi- 
static channel model as well as the standard technique to com- 
pute the FER performance of a scheme on that channel are 
briefly described in Section [TT] In Sections iHlWl we derive 
expressions based on which a threshold value, characteristic 
of the system, can be determined that allows the computation 
of a simple yet accurate approximation to the FER. Numerical 
results are presented in Section|Vl] whilst a method to produce 
performance evaluation plots is discussed in I VIII Finally, the 
main conclusions of our work are summarized in Section [Villi 



II. System Model 

If x is a frame of symbols transmitted over a quasi-static 
fading channel at a particular time instant and y is the receive 
frame, the input-output relationship of the channel is given by 



hx. 



(1) 



The instantaneous fading coefficient h is a zero-mean, circu- 
larly symmetric complex Gaussian random variable with vari- 
ance a 2 = 1, whilst n is a sequence of zero-mean, mutually 
independent, circularly symmetric complex Gaussian random 
variables with variance A/q. Note that h is constant for the du- 
ration of the transmit frame but changes independently from 
frame to frame. 

The quality of a quasi-static fading channel is character- 
ized by its corresponding average receive signal-to-noise ratio 
(SNR). In particular, if £ s is the energy per transmit symbol 
and 7= \h\ 2 £ s /Afo is the instantaneous receive SNR, the av- 
erage SNR, 7, at the input of the receiver is given by 

7 = E [7] = E [l^l 2 ] (£ S /M ) = a 2 (£ S /JV ) = £ S /N~ , (2) 

where E [.] denotes the expectation operator. 

The average FER on a quasi-static fading channel, denoted 
as P^(j), can be computed by integrating the FER in additive 
white Gaussian noise (AWGN), represented by P^(j), over 
the fading distribution [10] 



(3) 



Now, the fading magnitude \h\ has a Rayleigh distribution, 
so that the instantaneous value of 7 is chi-squared distributed 
with two degrees of freedom [10], i.e., 



Piil) = (1/7) e" 



-7/7 



for 7 > 0. 



(4) 



Although (0) is an exact expression for P~* (7), its eval- 
uation could prove difficult depending upon the transmission 
technique under consideration. In the following section we de- 
scribe a simple approach that is often used to bound the FER 
performance of a communication scheme. 



III. Approximation to the Frame Error Rate 

Given an arbitrary SNR threshold 7^, we can rewrite the ex- 
pression for the average FER on a quasi-static fading channel 
as follows 

P?{l) = P(error|7 < j w )P(j < j w ) 

+ P(error|7 > 7™)P(7 > J w ). 

The SNR threshold j w , which we refer to as the waterfall 
threshold, is used to divide the range of SNR values into a 
low-SNR region and a high-SNR region. A common approach 
to simplify (0 is to use a trivial bound of the form 



(5) 



P(error|7 < -f w ) < 1 



(6) 



for the low-SNR region and a conventional union bound for 
the high-SNR region. This approach has been used to bound 
the performance of many transmission schemes on quasi-static 
fading channels, including convolutional codes [4] and turbo 
codes [5]. Nevertheless, the value of "f w needs to be chosen 
appropriately to make the bound as tight as possible. The op- 
timization process presented in [4], [5] and [8] produced quite 
tight bounds but showed that there is room for improvement. 

Turbo codes have also been considered in [9] and [11]; the 
authors further simplified (0 by assuming that 

P(error|7 < 7^) w 1 and P(error|7 > -f w ) w 0. (7) 

Substituting (0 into (0 gives an approximation to the FER, 
denoted as P^ililw)' In particular, 

if (7) * P(7 < 7») 







(8) 



= 1 - e" 7 ™/^ 



It has been shown in [9] that P{? (7,7™) very accurately de- 
scribes the actual FER of a turbo code using a long interleaver 
on a quasi-static fading channel, if the waterfall threshold j w 
is set to be equal to the decoder convergence threshold jth- 
The convergence threshold of iterative schemes, such as turbo 
codes, can be determined using extrinsic information (EXIT) 
chart analysis [12]. 

Motivated by the work of Bouzekri and Miller [4], [5], El 
Gamal and Hammons [9] and Rodrigues et al. [1 1], we assume 
that for any transmission scheme, whose conditional FER can 
be reasonably described by there is a value of j w for 
which the average FER, Pe (7)> can be accurately approxi- 
mated by P,? (7,7to)- Based on that assumption, we derive an 
exact expression for the waterfall threshold in the following 
section. 

IV. Evaluation of the Waterfall Threshold 

Let £ denote the absolute difference between the actual 
frame error probability P®(pf) and the approximated frame 
error probability pQ(fy,^ w ), i.e., 



(9) 



where 7 can be any nonnegative real number. In the previous 
section we indicated that we expect the approximated FER of a 
transmission scheme on a quasi-static fading channel to closely 
represent the actual FER for a very wide range of 7 values, 
provided that an appropriate value for the waterfall threshold is 
chosen. Consequently, if P®(j,j w ) perfectly coincides with 
iv (t)> expression (0 simplifies to 



P?{l)~ P e Q (7,7«0 = 0. 



(10) 



In this section we derive an expression for the waterfall thresh- 
old j w under the assumption that e = 0, while in Section 
IvTl we compare our analytic approach to simulation results in 
order to test the validity of our assumption. 

We set A = 1/7 and express P®{^) and P® (7, 7^) as 
functions of A, i.e., P®(\) and pQ(\,j w ) respectively. The 
change of variable will not have any effect on s, hence 



= P«(A)-P«(A,7 W ) = 0, 



(11) 



for all values of A > 0. Equivalently, the area under the graph 
of P L ?(A) should be equal to the area under P®(\,j w ), for 
A e [0 ... A], where A — > 00. Consequently, we can write 

KmjjT P <?(A)dA-^ P?(\, lw )d\^ =0. (12) 

Using and (0]l, we expand the first integral in (fT2t into 

pA pA poo 

/ P®(\)d\ = / / P e G ( 7 )pA(7)^A 

JO JO JO 

Pe(l) / Ae" A7 dAd 7 . 



The frame error probability of a transmission scheme over an 
AWGN channel can also be expressed as P e G (7) = 1 — P^il), 
where Pfil) is the probability of successful frame detection. 
Consequently, we can rewrite < TT~3T > as 



P?{\)d\ 



Xe-^'d\d 1 - 

(14) 

Pdh) / Xe^dXdj. 
Jo 



Taking into account that [13] 



A 
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Ae~ A7 dA = — 
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,-A-y 



(1 + A 7 ) 



and 



\ e - x ~<d\dj = A, 

10 Jo 

the first integral in ( fl2b assumes the form 



(15) 



(16) 



P e Q (A)dA =A- 
Jo 
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,-A7 



(17) 



+ / P d(l)— ^(1 + A 7 )d7- 



The second integral in ( fTSl i can be evaluated as follows 



Jo Jo 



-A 



1-e 



rfA 



1 



(18) 



7™ Tt« 

If we substitute dTTb and ( fT~8T > into ( fT2b . we observe that 
terms A and —A cancel each other out. Furthermore, if we take 
the limit as A — > 00, all terms containing e~ A are eliminated 



since e 



0. The remaining terms give 



I), 



which is equivalent to 



(19) 



(20) 



We have thus shown that, under assumption (flOl l. the waterfall 
threshold is inversely proportional to the area under a curve, 
which is defined by the probability of successful frame detec- 
tion in AWGN normalized by the squared instantaneous SNR. 

Depending on the expression for the detection probability 
Pdij), a closed-form solution for j w for a particular transmis- 
sion technique may not exist. In that case, 7™ can be evaluated 
either via numerical integration if there is an exact representa- 
tion of P<p (7) or via Monte Carlo simulation if Pj* (7) cannot 
be accurately evaluated. In the following section we revisit 
(f20b to obtain a more practical expression for j w when Monte 
Carlo simulation is required. 

V. Practical Computation of the 
Waterfall Threshold 

Let y be the actual SNR value for which Pf (7) = if 
7 < 7' but P2{l) > otherwise. Based on the definition of 



we can rewrite d20b as follows 



(7) 



el', 



or, equivalently, 



7« 



1 

7 



p?(.i) 



7' 



r 



(21) 



(22) 



It becomes evident in Fig. [T] that as 7 grows, function 
•P<P (7V7 2 gradually approaches I/7 2 , which slowly con- 
verges towards zero. The advantage of (l22l over (l20b is that 
P^{l)ll 2 converges to zero much faster than Pff ("f) / "f 2 , 
hence an accurate value for j w can be obtained by considering 
only a limited low-SNR range of integration. 

Let us now consider the case when Monte Carlo simulation 
is used to measure the FER performance of a transmission 
scheme in AWGN. We assume that the SNR values 7$, with 
i = 1, 2, . . . , N, are equally spaced and ordered, while the FER 
is -P e G (7i) = 1 for i < k but P^ili) < 1 otherwise. Elaborating 




Fig. 1, Normalized probabilities in AWGN. In this example, we have con- 
sidered an input frame length of 512 bits and a system using a recursive 
convolutional code with generator polynomials (1,17/15) in octal form. 



on d22l . we can obtain the following equivalent expression for 
discrete SNR values 




(23) 



Substituting ( 1231 into dD gives us an approximation of the av- 
erage FER of the transmission scheme on a quasi-static fading 
channel. 

Note that the complexity of the threshold-based FER com- 
putation using (0 and either ( f20l > or (|23l is markedly less than 
that of the exact FER computation based on (01, as we will 
now demonstrate. Let us consider the case when the frame er- 
ror probability in AWGN, P^ (7), is known for A~ values of 7 
and our objective is to compute the FER for a quasi-static fad- 
ing channel, pQfi), for M values of 7. The threshold-based 
approach involves the evaluation of the waterfall threshold, 
which has a computational complexity proportional to N, fol- 
lowed by the calculation of the FER approximation, which has 
a computational complexity proportional to M. Consequently, 
the overall complexity of the threshold-based FER computa- 
tion is of order 0{N + M). In contrast, computation of the 
exact FER requires A^ multiplications for each value of 7, 
resulting in a complexity order of 0(NM). 

VI. Numerical Results 

In this section we compare analytical to simulation results 
for transmission over quasi-static fading channels. We consider 
both uncoded and coded binary phase shift keying (BPSK); 
coded transmission uses either a rate 1/2 recursive systematic 
convolutional code with octal generator polynomials (1,17/15) 
or a rate 1/3 turbo code with generator polynomials (1,5/7,5/7). 
The input frame length L is either 256 bits or 1024 bits. 
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Fig. 2. Frame error rate performance of various transmission schemes on a Fig. 3. Performance plots for the convolutional and turbo codes under con- 
quasi-static fading channel for input frame lengths of 256 and 1024 bits. sideration. 



The waterfall threshold of uncoded BPSK was computed 
numerically using d20l ). where the probability of successful 
frame detection is captured by [14] 



(24) 



Here, Q(x) is the tail integral of a standard Gaussian density 
with zero mean and unit variance, defined as 



Q(x) 



(v 



(25) 



The threshold j w was found to be 5.782 dB for L = 256 and 
7.083 dB for £ = 1024. 

In the case of coded transmission, the probability of erro- 
neous frame detection in AWGN was obtained running Monte 
Carlo simulations for a limited range of low SNR values (e.g., 
7 G (0, 1] for turbo coding) and a small number of channel 
realizations (a few thousands at most). The waterfall threshold 
was then calculated using j23l . In particular, when convolu- 
tional coding is employed it was found that ~f w = —0.983 
dB for L = 256, while <y w = 0.023 dB for L = 1024. When 
turbo coding is used, however, the frame length appears to 
have minimal impact on the waterfall threshold; its value was 
found to be j w = -4.401 dB when L = 256 and j w = -4.312 
dB when L = 1024. Hence, we expect that the input frame 
length will not significantly affect the FER performance of 
the turbo code. 

Once the waterfall threshold of a scheme has been com- 
puted, we substitute it into ^ to obtain an analytical expres- 
sion for the approximated average FER for the quasi-static fad- 
ing channel. The curves of the approximated FER expressions 
for the systems under investigation are compared to simula- 
tion results in Fig. [2] Observe that the analytic technique very 
closely approximates (within 0.4 dB) the simulation results in 
the various scenarios. Hence, both coded and uncoded, itera- 
tive and non-iterative systems can indeed be characterized by 
waterfall thresholds based on which tight approximations for 



their frame error probability can be derived. It is also inter- 
esting to note that, as expected, the FER performance of the 
turbo code remains unaffected by the input frame length, or 
equivalently the interleaver size. The same behavior has also 
been reported in [5], [11]. 

Based on the comparison between analytical and simulation 
results, we conclude that the technique presented in this paper 
accurately estimates the FER performance of BPSK transmis- 
sion schemes over quasi-static fading channels. Similar results 
can be easily obtained using the same approach for higher 
order modulations. 

VII. Discussion on Performance Evaluation 

An insight into the relative performance of two or more 
transmission schemes on quasi-static fading channels can also 
be obtained by plotting their normalized probabilities of suc- 
cessful frame detection in AWGN, i.e, P2{l)ll 2 ^ an d com- 
paring the areas under the corresponding curves. In particular, 
the larger an area is, the smaller is j w and, consequently, the 
better the FER of the scheme under consideration is expected 
to be, according to <[8j and S20) , 

The normalized probabilities P2{l)ll 2 of the coded BPSK 
schemes that we considered in the previous section, namely 
the rate 1/2 convolutional code and the rate 1/3 turbo code, 
have been plotted in Fig. [3] We observe that as the input frame 
length of the convolutional code increases, the area under the 
normalized probability graph reduces and, as we have already 
seen in Fig. [21 the FER performance of the convolutionally 
coded scheme degrades. As anticipated, the area under the 
normalized probability curve of the turbo code is clearly larger 
than that of the convolutional code, hence the former scheme 
yields better FER performance on quasi-static fading channels. 
Most importantly however note that, in the case of the turbo 
code, as the length of the input frame increases, the shape of 
the curve changes such that its peak shifts leftward but at the 
same time, moves upward. For L large, it is expected that the 
peak will reach its maximum value, which is a point on I/7 2 ; 




y 



Fig. 4. Performance plots of a turbo code with generator polynomials 
(1,5/7,5/7), for various interleaver sizes. 

notice that the curve 1 /j 2 corresponds to the ideal case where 
P9{p/) = 1 for all values of 7. 

This trend is more clearly illustrated in Fig. |U where the 
normalized probability of the turbo code has been plotted for 
various input frame lengths or, equivalently, interleaver sizes. 
In all cases, the exact log-MAP algorithm was used, while all 
probabilities were measured after 8 decoding iterations. We 
can infer from Fig. [4] that the probability of successful frame 
detection of the turbo code becomes a step function when 
very large interleavers are used, i.e., L — ► 00. In particular, 
P d{l) = for 7 < -fth whereas P^il) = 1 for 7 > j th ; 
note that -f t h corresponds to the convergence threshold of the 
iterative decoder [9]. Using (|20| |. we can obtain the waterfall 
threshold of the turbo code for L — > 00, as follows 

lw=(^J ~2 d ^ =lth- (26) 

Therefore, our approach is in complete agreement with the 
findings of El Gamal and Hammons [9], i.e., the approximated 
FER expression is an accurate representation of the actual FER 
of a turbo code using a long interleaver on a quasi-static fading 
channel, if the waterfall threshold is set to be equal to the 
convergence threshold of the iterative decoder. 

VIII. Conclusions 

In this paper, we have considered various transmission 
schemes, both uncoded and coded, iterative or non-iterative, 
over quasi-static fading channels and we have demonstrated 
that a waterfall threshold can be used to characterize them. We 
have provided an exact interpretation of the waterfall thresh- 
old, based on which an accurate approximation of the frame 
error rate can be obtained. Finally, we have analytically con- 
firmed that our approach is in agreement with the literature, 
when turbo codes are considered; in particular, we have shown 
that the waterfall threshold indeed coincides with the conver- 
gence threshold of the corresponding iterative decoder, when 
long interleavers are used. 
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